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INTRODUCTION 


An  often  studied  and  important  series  of  rare  earth-transition 
metal  compounds  are  those  that  form  in  the  cubic  Laves  or  Cl 5 phase. 
In  this  structure  there  are  two  rare  earth  and  four  transition  metal 
atoms  in  the  unit  cell.  The  rare  earth  atoms  have  tetrahedral  (T^) 
site  symmetry,  i.e.,  the  positions  of  the  carbon  atoms  in  the  diamond 
structure,  while  the  lighter  metal  atoms  have  03,3  symmetry.  (See 
Figure  1.)  In  the  compounds  with  Fe,  the  rare  earth  spins  are  taken 
to  be  parallel  to  one  another  and  antiparallel  to  the  iron  spins,  a 
picture  which  is  consistent  with  all  revelant  experiments. Rare 
earth  compounds  with  iron  in  the  Laves  (C15)  phase  are  strongly 
magnetic  to  well  above  room  temperature,  and  TbFe2  has  the  largest 
known  magnetostriction  and  magnetic  anisotropy  at  room 
temperature. (2)  These  last  two  properties  derive  from  the  rare 
earth,  and  accordingly  disappear  for  those  magnetic  compounds  like 
YFe2  and  GdFe2 , where  the  rare  earth  atom  has  no  orbital  moment.  It 
has  also  become  increasingly  clear (3)  that  the  magnetostriction  for 
these  compounds  is  very  anisotropic.  For  example,  TbFe2  which  has 
its  easy  magnetization  direction  parallel  [111]  develops  strains  of 
2x10”3,  while  DyFe2 , where  [100]  is  easy  exhibits  strains  of  10”5  or 
less.  Since  rare  earth  environments  in  these  compounds  are  identical 
in  every  way,  this  anisotropy  is  expected  to  be  a common  feature  of 
all  the  compounds.  Thus  we  take  as  a fact  that  ^m>>Aion  for  Laves 
phase  rare  earth  intermetallics . This  is  unusual  for  cubic  metals. 
The  cubic  ferromagnets  Fe  and  Ni  for  instance,  show  no  such  effect. 

We  expect  little  anisotropy  in  the  elastic  constants.  If  the 
magnetoelastic  coupling  energies  Bi  (100  easy)  and  B2  (111  easy)  are 
the  same,  the  predicted  (4)  ratio  of  strain  in  TbFe2  to  that  in  DyFe2 
is  about  two;  the  experiments  give  a ratio  more  like  one  hundred. 

^Taylor,  K.  N.  R. , Advances  in  Physics,  20,  551  (1971).  This 
article  gives  an  extensive  review  of  pertinent  experiments 
relating  to  rare  earth  intermetallics. 

(2)ciark,  A.  E.  and  Belson,  H.  S.,  17th  Conference  on  Magnetism  and 
Magnetic  Materials,  Chicago,  Illinois,  AIP  Conference  Proceedings 
No.  5,  1498  (1972);  (Phys.  Rev.  B 5,  3642  (1972))  Clark,  A., 
Cullen,  J. , and  Sato,  K. , Proc.  20th  Conf.  on  Magnetism  and 
Magnetic  Materials,  AIP  Conf.  Proc.  24,  670  (1975). 

Clark,  A.  E.,  19th  Conference  on  Magnetism  and  Magnetic  Materials 
Nov  1973,  AIP  Conference  Proc.  18,  1015  (1974) 

(4^The  prediction  is  based  on  the  use  of  equations  (9)  and  (12)  of 
this  paper. 
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In  this  paper,  we  describe  a theory  for  the  spontaneous 
distortions  of  the  atoms  in  the  cubic  Laves  phase  lattice  due  to 
coupling  of  the  crystalline  electrostatic  fields  with  the  total 
angular  momentum  of  the  4f  electrons. '^)  in  the  usual  picture  of 
magnetostriction,  this  distortion  is  a uniform  strain  in  either  of 
the  two  external  strain  modes:  2eZz  " exx  “eyy»  exx  “ eyv'  a 

distortion  along  100  (*ioo)'  and  exy»  z\z>  exz»  a distortion  along 
111  (Am)  where  x,  y,  z refer  to  a cartesian  set  of  axes  parallel  to 
the  [100]  directions.  We  show  here  that  in  the  C15  structure,  where 
there  are  two  rare  earth  atoms  in  the  unit  cell,  the  crystalline 
field-angular  momentum  coupling  gives  rise  to  an  internal  distortion 
which  in  turn  is  coupled  to  an  external  strain. The  results  are  a 
predicted  internal  distortion  of  the  rare  earth  sublattice  and  an 
additional,  extra-ordinary  contribution  to  the  external  strain.  This 
additional  strain  offers  a natural  explanation  for  both  the  size  and 
anisotropy  of  the  magnetostriction  in  the  Laves  phase  intermetallics . 

In  the  succeeding  sections  of  this  paper  we  give  the  physical 
picture  of  the  process  by  which  a C15  crystal  distorts  in  the 
presence  of  a net  magnetic  moment.  We  then  calculate  the  magnitudes 
of  the  distortions  as  a function  of  crystal  field  energies,  internal 
and  external  moduli,  and  the  shape  of  the  4f  electron  charge 
distribution.  We  give  numerical  estimates  of  the  sizes  of  the 
predicted  effects  based  on  a point-charge  model  and  experimental 
values  for  the  moduli. 

1.  The  Physical  Picture 

The  way  in  which  a spontaneous  distortion  occurs  can  be  seen  as 
follows:  Suppose  the  magnetic  moment  is  parallel  to  a [111] 

direction.  Because  of  strong  spin-orbit  coupling  the  charge  density 
of  every  rare  earth  atom  takes  on  a very  special  shape  for  4f 
electrons.  For  Tb,  for  instance,  this  density  is  pancake-like, 
(oblate)  lying  in  a plane  perpendicular  to  the  [111]  direction,  as 
shown  in  Figure  2.  The  charge  cloud  of  the  atom  labelled  A in  this 
figure  comes  closer  to  the  three  near  neighbor  atoms  labelled  B'  than 
to  its  near  neighbor  at  B.  Because  of  a net  coulomb  attraction  of 
this  charge  with  the  positive  cores  of  its  neighbors,  the  A atom 
moves  toward  the  B'  atoms,  i.e.,  down  in  the  Figure.  Exactly  the 

(5) The  Phenomenology  of  Magnetostriction  was  first  presented  in 
R.  Becker  and  W.  Doring,  Ferromagnetismus  (Verlag  Julius 
Springer,  Berlin,  (1939)  P.  132,  145.  See  also  C.  Kittel, 
"Introduction  to  Solid  State  Physics"  2nd  edition  John  Wiley  and 
Sons,  Inc.  (New  York).  Reference  9 of  this  paper  describes  a 
microscopic  calculation  of  magnetostriction  of  rare  earth  metals. 

; The  possibility  of  internal  distortions  coupled  to  external  strain 
is  mentioned  in  E.  R.  Callen  and  H.  B.  Callen,  Phys.  Rev.  129,  578 
(1963). 
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opposite  occurs  at  B;  this  atom  moves  up  in  response  to  the  net 
attraction  of  A'  atoms  over  the  A atom.  Figure  3 shows  the  situation 
when  the  magnetization  is  along  [001] , the  easy  axis  for  DyFe2  and 
HoFe2 . The  charge  cloud  remains  equidistant  from  all  its  neighbors, 
and  no  distortion  can  occur.  We  now  describe  the  calculation  of  this 
distortion  and  the  consequent  external  strain. 

2.  Calculation  of  the  Strains 

The  Electron-Lattice  Interaction  - Consider  a 4f  electron 
associated  with  a rare  earth  atom  at  the  origin.  Its  wave  function 
is  changed  from  its  atomic  form  by  the  crystalline  potential  of  the 
nearby  rare  earth  and  iron  atoms.  The  electron  is  described  by  its 
total  angular  momentum  vector  J.  An  effective  Hamiltonian  for  a 
given  site  can  be  derived  for  the  crystalline  field  by  replacing  the 
electron  coordinates  by  their  angular  momentum  operator 
equivalents. (?)  The  first  step  is  to  obtain  the  correct  form  of  the 
crystalline  field. 

For  a site  with  tetrahedral  symmetry, 

OOO  AAA 

Vcf  = u2  (x  +y  +z  ) + u3  xyz  + u4  (x*+y  +z^)  + ... 

x,y,z  are  the  dimensionless  electron  coordinates  normalized  by  the 

lattice  parameter,  cIq.  Now  because  of  parity,  the  average  of  the 

second  term  of  Vcf  vanishes;  i.e.,  there  is  no  effect  from  the 

odd-order  term.  Suppose,  however,  that  the  ion  moves  to  a new 

position  with  coordinates  UxA>  uyA*  uzA  say,  with  no  change  in 

size  or  shape  of  the  unit  cell.  The  energy  of  the  4f  electron  then 

depends  on  u.  To  calculate  Vcf  (u)  we  proceed  as  follows:  In  the 

electron-ion  potential  the  electron  coordinates  x,y,z  are  replaced  by 

x-ux,  y-Uy,  z-u2.  We  must  also  take  into  account  the  movement  of  the 

neighboring  rare  earth  sublattice.  Thus  the  ion  coordinates  Xn,Yn,Zn 

are  replaced  by  Xn-uxn,Yn-un,  Z-uzn.  At  this  stage  we  make  the 

assumption  that  u^  = ub  for  all  near  neighbors  to  an  A site. 

Essentially,  the  reason  this  simplification  works  is  that  the  near 

neighbors  of  a given  rare  earth  on  one  sublattice  (A)  are  all  members 

of  the  other  sublattice  (B) . Since  the  electron-ion  potential 

contains  x and  xn  only  in  the  form  x-xn,  the  distortion  coordinates 

appears  only  as  the  difference  in  distortions  on  A and  B sites. 

Define  these  coordinates  to  be  u ,u  ,u  . We  therefore  have 

x r y z 

Vcf  = v2A  £(x-ux)2  + (y-uy)2  + (z-uz)2l  + v3A  (x-ux)  (y-Uy)  (z-uz) 


7TJ 

Stevens,  K.  W.  H. , Proc.  Phys.  Soc.  (London)  65,  209  (1952) 
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plus  terms  of  higher  order  in  x-ux.  The  coefficients  V2,V3,  etc., 
are  lattice  sums;  in  view  of  the  special  importance  of  the  third 
* order  term  it  may  be  worthwhile  giving  the  expression  for  V3  on  the 

point-charge  model; 

..  A _ lc„_2  VXnYnZn  3 . , „ 

v3  = -15Ze  7 — aQ  1(a) 

n R 
n 

Here  Z is  the  effective  number  of  charges  on  a rare  earth  ion  and  e 
is  the  electron  charge.  The  sum  is  over  all  near  neighbors.  Now 
consider  atom  B at  (1/4,  1/4,  1/4).  The  form  of  the  crystal  field  is 
the  same,  but  u appears  with  opposite  sign 

V^f  = v2B  |u+ux)2  + (y+uy)2  + (z+uz)2J+  v3B  (x+ux)  (y+uy)  (zHiz)  + . 

The  different  environments  on  A and  B atoms  have  no  effect  on  V2;  V2A 
= v2.  For  the  third-order  term  however,  V3A  = -V3B  as  can  be  seen 
by  examining  the  point  charge  expression  for  V3.  For  an  A atom  at 
the  origin,  the  sum  is  over  neighboring  rare  earths  at  1/4,  1/4,  1/4; 
1/4,  -1/4,  -1/4;  -1/4,  -1/4,  1/4;  -1/4,  1/4,  -1/4,  in  units  of  the 
lattice  parameter.  For  a B atom  (considered  similarly  located)  , the 
corresponding  A neighbors  would  be  at  -1/4,  -1/4,  -1/4;  -1/4,  1/4, 
1/4;  1/4,  -1/4,  1/4;  1/4,  1/4,  -1/4. 

So  far,  we  have  only  considered  a possible  energy  lowering  by 
internal  distortion,  u.  We  can  also  expect  that  the  symmetry  change 
due  to  this  distortion  may  drive  an  external  strain  as  well,  i.e.  , an 
extraordinary  magnetostriction  may  appear,  causing  a shape  change  of 
the  unit  cell.  We  can  take  account  of  external  strain  and  thereby 
account  for  both  ordinary  and  extraordinary  magnetostriction  by 
writing: 


« 


Xn  " Xn°  (1+exx> 


+ Yn°  exy/2 


Z ° e /2 
n xz' 


etc. , for  the  coordinates  of  the  atom  at  the  nth  site.  Because  this 
strain  lowers  the  symmetry  at  A and  B sites,  there  appear  lower-order 
terms  in  the  crystal  field.  The  field  at  an  A site  for  example  is 
now : 


Vcf  (u'e)  = vl"  (x  ,eyz+y,eYZ+z,exy)  + v3A  x'y’z' 


xy 


+ v0A  (z,2+...)  + v"  , (x'y'e  + ...) 


'2,1 


5 


1(b) 
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!• 


Here,  as  in  the  constrained  case,  x'  = x-ux,  etc.  On  the  point 
charge  model 


f z 


n 


Rn 


5“  ao 


(2) 


>3a  is  given  by  Equation  (la)  . 


The  sum  is  over  the  unstrained  coordinates  of  near  neighbors  to  an  A 
site.  It  is  clear  that  in  general 


I 


i: 


viA  = v*,  i = 2n 

vj.A  = -v±B,  i = 2n  + 1 (3) 


The  electron-lattice  energy  density,  or  the  magneto-elastic 
Hamiltonian  if  operator  equivalents  are  used,  is  the  sum  of 
contributions  from  each  sublattice  plus  the  elastic  energy  density: 


E + E 
e me 


<VAC>  + 
cf 


<V®->  + 7 C.  (e  2+e  2+e  2) 

cf  2 44  xy  xz  yz 


+ 2 (Cll"C12^  (evv+eww+eiL)  + 7 k (u^+u,2,+u?) 


xx  yy  zz 


I (4) 

The  C's  are  the  usual  elastic  moduli,  k the  appropriate  internal 
modulus.  Volume  strains  do  not  concern  us,  and  have  been  dropped. 
Also,  we  have  assumed  isotropic  internal  forces.  There  are,  of 
course,  many  more  internal  lattice  modes,  but  they  are  not  coupled  to 
uk  an<*  £i j / and  so  need  not  be  added  to  the  energy  expression.  It  is 
more  convenient  to  express  the  resulting  equations  in  terms  of  the 
eij,  uk  and  crystal  field  energy  densities  Ei#  where  Ei  = p Vj.  The 
final  expression,  then  for  the  magneto— elastic  energy  density, 
suppressing  terms  which  do  not  contain  uj^  or  or  which  give  zero 
in  the  4f  ground  state,  J 


E 


me 


2E1  (uxeyz+*“)  + 2E3  (ux  <YZ>  + •••) 

+ 2 E,  1 (<X2>  e +...)+  2 E,  , (<XY>  e 
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Actually  the  term  in  Ei  in  this  expression  should  more  properly  be 
included  below  in  the  elastic  energy,  since  it  does  not  depend  on  the 
angular  moments  of  the  4f  electrons;  the  coupling  it  represents 
persists  above  the  Curie  point  and  should  exist  even  in  YFe2,  GdFe2 
where  there  are  no  orbital  moments.  It  is  in  fact  a property  of  all 
lattices  of  the  diamond  or  zinc-blend  structure  that  the  internal  and 
external  strains  are  so  coupled.  Thus,  a proper  calculation  of 
presents  the  same  order  of  difficulty  as  the  calculations  (not 
attempted  here)  of  the  elastic  constants  and  internal  moduli  defined 
below. 


Finally, 


E = i cAA  (e  2+e  2+e  2)  + i (C,,-C,,J  (e  2+e  2+e  2) 
e 2 44  xy  xz  yz'  2 v 11  12'  ' xx  yy  zz' 

+ 7 k ^ux2+uy2+uz2^  ' 


The  total  strain-dependent  energy  is  therefore  E Ee  + Eme.  The  E's 
all  have  units  of  energy,  k is  the  internal  modulus.  We  have 
written  E^  = E^A,  etc.,  to  simplify  notation. 

The  Equilibrium  Strains  - The  result  of  minimizing  E with 
respect  to  the  u's  and  e*s,  i.e. , 


kux  - 2E^  e = 2E3  <YZ>,  etc. 


C44  Eyz  ' 2E1  “x  ’ 2E2,2  <VZ>'  etc' 


(Cll‘C12)  exx  2E2,1  <X  >*  etc‘ 
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Note  that  ux  is  coupled  to  eyz,  uy  to  ezx  and  uz  to  exv 
not  coupled  to  internal  modes . * 


The  are 


We  deal  with  more-than-one  electron  configurations  in  the  standard 
way:  replace  xy  for  example  by  E Xj^,  where  the  sum  is  over  all  4f 

electrons,  and  then  using  the  operator-equivalent technique. 

Hence : 


<XY> 


a 


(JxJy+JyJx) 


(ID 


a's  for  various  rare  earths  are  tabulated. (5)  <r2f>  is  the  average 

of  the  square  of  the  4f  electron  radius,  which  is  also  tabulated. ^ 

At  temperatures  well  below  the  Curie  point  this  expression 
(Equation  11)  may  be  replaced  by  its  zero-temperature  valued' 


J J + J J 
x y y x 


(j2-j) 


fxTy  / 


where  y's  are  the  direction  cosines  of  the  magnetization. 
Thus 


t 

t 


n 


! 


i 


I 

i 


<XY>  = J(J-l/2)a<r2f>  . (12) 
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' 1 


For  the  lattice  to  be  stable  under  these  strains  kC.. 

44 


> 4E, 


3.  Numerical  Estimates 


We  now  proceed  to  a discussion  of  the  sizes  of  the  predicted 
effects.  Although  the  strain  modulus  is  known  for  Laves  phase 
compounds , (10)  the  internal  modulus  is  not.  In  principle  it  could  be 
obtained  from  measurements  of  the  energy  at  long  wavelengths  of  a 
special  "optic"  mode  in  the  paramagnetic  phase  of  these  compounds. 

The  special  mode  we  mean  is  that  corresponding  to  the  two 
inequivalent  rare  earths  oscillating  180°  out  of  phase  with  one 
another  while  the  iron  sublattices  sit  still.  It  can  be  shown  (See 
Appendix)  rather  simply  that  this  sort  of  oscillation  is  in  fact  a 
normal  mode  in  the  cubic  Laves  phase  if  all  forces  are  assumed 
isotropic.  Indeed,  this  is  the  reason  we  could  neglect  any  possible 
iron-atom  displacement  in  the  preceding  discussion;  such 
displacements  just  do  not  couple  to  the  internal  mode  we  are 
discussing.  As  we  show  in  the  appendix  the  modulus  k can  be  written 
in  terms  of  the  iroh-rare  earth  and  rare  earth-rare  earth  force 
constants : 


k 


2KR,R  + 


4K 


F , R 


Although  the  point  charge  model  for  calculating  the  crystal 
field  parameters  and  their  strain  derivatives  is  not  expected  to  be 
accurate,  we  use  it  as  a means  of  estimating  ratios  of  these 
parameters.  That  is,  we  calculate  E3  and  E^  in  terms  of  E2r2  and 
determine  E2,2  from  magnetostriction  values  of  a metallic  ferromagnet 
with  one  atom  per  unit  cell.  On  the  point  charge  model,  E3  = 4OE2  2 
and  E^  = — (3/4) E^  2*  Assuming  4E^2/kc  <1,  e reduces  to 


' I 

:1 

1 J 


I 

( 


e 


yz 


(60E_  _/k+l)  <YZ>  . 

C44 


To  obtain  an  idea  of  the  size  of  the  internal  modulus  required  to 
make  the  extra-ordinary  term  dominant  in  this  expression,  we  set  the 
ordinary  part  equal  to  a value  typical  of  cubic  ferromagnets : 


2E 


(e  ) , . 

yz  ordinary 


<YZ>  £ 10"5 


'44 


(10)  Clark,  A.  E.  and  Savage  H.  T. , IEEE  Trans.  Sonics  and 
Ultrasonics  SU-22,  50,  January  1975. 
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From  tables (8)  of  Steven's  factors  we  find  that,  for  Tb^)  at  low 
temperatures , 


<YZ>  = .33  x 10  J . 

2 ° 2 o 

Here  we  have  used  r f = .185(A)  and  aQ  = 8 A.  (The  radius  was 
obtained  from  reference  (6)  and  the  lattice  constant  from  reference 
(1).)  Finally,  using 


C.A  = 2x10  ergs/cm' 


we  get 


9 3 

E2  2 = 10  ergs/cm 

Thus,  for  the  "extra-ordinary"  part  of  the  strain  to  dominate, 
and  give  the  experimental  result  for  TbFe2,  we  require  E2  2/k~l> 
which  means  the  normal  frequency  of  the  internal  mode  in  the 
paramagnetic  phase  would  have  to  be  lO*2  cps  or  lower.  Such 
frequencies  are  typical  of  optic  modes  in  metals  with  more  than  one 
atom  per  cell.  In  sum,  the  size  and  anisotropy  of  the  magneto- 
striction in  Laves  phase  magnets  can  be  accounted  for  on  the  basis  of 
the  theory  given  here,  culminating  in  Equation  (13) . The 
magnetostriction  is  predicted  to  be  positive  for  TbFe2,  (a  is 
negative)  on  the  basis  of  Equation  (12)  if  the  extraordinary  term 
dominates.  This  accords  with  experiment. (2) 

The  prediction  of  an  internal  distortion  can  be  checked  by  X-ray 
measurements,  though  the  analysis  is  complicated  by  the  distortion 
caused  by  the  external  strain.  We  want  to  make  a remark  about  the 
direction  of  the  internal  distortion.  It  follows  from  Equation  (14) 
that  the  distortion  is  parallel  to  the  magnetization  when  the  latter 
lies  along  a [111]  axis;  it  vanishes  when  the  moment  is  parallel 
[100] . When  the  moment  is  in  the  xy  plane,  the  A and  B sublattices 
are  alternately  distorted  positively  and  negatively  along  the  Z axis. 
For  a general  direction  of  magnetization  the  two  sublattices  undergo 
antiparallel  distortions  but  the  axis  of  the  internal  distortion  is 
not  along  the  moment  direction  whereas  external  distortion  is  along 
magnetization  direction.  This  strange  behavior  is  consistent  with 


( 8 )  Freeman,  A.  J.  and  Watson,  R.  E.,  Phys.  Rev.  127,  2058  (1962). 

(9)  Tsuya,  N.  Clark,  A.  E. , and  Bozorth,  R. , Proc.  International 
Conf . on  Magnetism,  Nottingham,  p.  250  (1964) . 
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the  physical  picture  we  gave  in  Section  1:  The  4f  charge  cloud 
always  tries  to  move  so  as  to  get  closer  to  its  nearest  set  of 
neighbors  (See  Figure  2) . 

Finally,  we  note  that  these  effects  should  also  occur  in  the 
isostructural  Co  and  Ni  C15  compounds  with  the  rare  earths  at 
temperatures  below  their  respective  transition  temperatures.  A 
similar  internal  distortion  will  also  take  place  in  the  magnetic 
phases  of  the  similar  C14  compounds,  with  the  added  complication  that 
the  latter  are  hexagonal,  so  that  there  are  more  than  two  external 
strain  modes. 


In  our  description  of  the  atomic  displacements  we  have  ignored 
distortions  of  the  transition  metal  atoms  and  taken  the  rare  earths 
to  move  equal  amounts  in  opposite  directions.  For  this  theory  to  be 
a complete  account  of  the  magnetostriction  requires  that  this  rare 
earth  movement  be  a normal  mode  of  vibration  in  the  absence  of  strain 
and  that  there  be  no  other  modes,  i.e.,  iron  sublattice  motion,  which 
couple  to  the  external  strain.  The  first  of  these  requirements  is 
satisfied,  as  we  show  in  the  Appendix.  The  second  is  not,  a 
situation  we  will  deal  with  in  a subsequent  publication.  The  rare 
earth  displacement  is,  we  feel,  the  most  important  one  for 
considerations  of  magnetostriction. 
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APPENDIX  A 

Our  purpose  here  is  to  show  that  we  are  justified  in  assuming 
that  n A i.e.,  that  is  no  coupling  to  distortions  of  the 

transition  metal  sublattices.  (That  all  the  rare  earth  atoms  of  a 
given  sublattice  undergo  the  same  distortion  follows  from 
translational  invariance.)  We  do  this  by  calculating  the  normal 
modes  of  the  cubic  Laves  phase  assuming  near-neighbor  coupling.  For 
the  set  of  solutions  which  vary  in  time  like  exp  (-iut)  in  the  limit 
of  long  wavelength  the  equations  of  motion  are 


-io2M  u^  =J]  J]  (u*  -U^  ). 

a a j 67a  a'6  6 a 


Where  a and  8 are  sublattice  indices  running  from  one  to  six,  and  i 
and  j refer  to  x,  y,  or  z components. 

The  force  constants  are  defined  as  : 

a8 


9 


where  the  sum  is  over  all  near  neighbors  to  an  a site  which  are  on 
the  8 sublattice.  Label  the  rare  earth  atoms  by  a.=  5 and  6.  If  V 
depends  only  on  the  distance  between  sites  then  K1^  = 0 if  i / j 


and  a = 5,  8=6. 

2 


Also  „ = K1J 


'5,8 


6,8 


a8 

for  8=1  through  4. 


Then 


-w"M  u~ p.  = (uiR-u1R)  K^a  (uj-uj) 


'5,8 


j B=l, 

4 


5,6 


8 


-“2"  u\  " K6^6 


j 6=1 


5,6 


Here  M is  the  rare  earth. mass.  .It  can  also  be  shown  from  the 
definition  of  that  so  that,  by  subtracting  the 

last  two  equations  5,8  6,6 
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-A  = 2K5(6  (ui6-uis)  - (u^-u1,)  E K5i6  . 

4 

Defining  K_  , = K and  1>K_  . = 4 K_  we  see  that  there 

DfO  K^K  3=1  r fR 

exists  a normal  mode  in  which  the  two  rare  earths  beat  against  one 
another  with  a frequency 


Further,  the  normal  coordinates  for  this  mode  have  the  property  that 
uA  = -uB,  ug  = 0 for  all  four  transition  metal  atoms.  This  is 
exactly  the  form  of  static  distortion  we  have  assumed  in  Section  2. 


I 


I 
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Figure  Captions 


Figure  1 

A drawing  of  the  atomic  arrangements  in  the  C15  structure.  The  R 
atoms  are  arranged  identically  to  the  carbon  atoms  in  the  diamond 
latti-ce*  Each  unit  cell  has  two  R (rare  earth)  and  four  T 
(transition  metal)  sites. 

Figure  2 

The  distortion  along  [111]  which  results  from  the  different  net 
Coulomb  attractions  of  the  two  rare  earth  ions  in  the  C15  unit  cell, 


Figure  3 

When  [100]  is  easy,  there  is  no  distortion;  the  charge  cloud  on  the 
rare  earth  remains  equidistant  from  all  its  neighbors. 
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